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We consider the pion matrix elements of the isoscalar and isovector combinations of the vector 
and tensor twist-two operators that determine the moments of the various pion generalised parton 
distributions. Our analysis is performed using partially-quenched chiral perturbation theory. We 
work in the SU(2) and SU(4|2) theories and present our results at infinite volume and also at finite 

\Q • volume where some subtleties arise. These results are useful for extrapolations of lattice calculations 

f^*) ' of these matrix elements at small momentum transfer to the physical regime. 
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Q ', Generalised parton distributions (GPDs) [H, H, H, 0] provide a uniquely detailed view of the structure of hadrons, 
I ■ unifying the information encoded in form-factors and parton distributions, and supplementing both. Ongoing ex- 
\ periments at DESY [H, 0, 0] and Jefferson Lab [1] (see Refs. [l(| for recent reviews) seek to learn about these 
fundamental quantities in deeply-virtual Compton scattering (DVCS) and related processes. Aspects of GPDs are 
also being investigated in QCD and phenomenological models. Since GPDs encode long distance hadronic structure, 
QCD analyses of them are necessarily based on non-perturbative methods such as lattice QCD. These examinations 
are complementary to the experimental efforts in that different facets of the GPDs can be accessed. Most experimen- 
tal efforts are focused on the proton and either measure the GPDs at constrained kinematics (through polarisation 
observables) or measure integrals of the GPDs (through DVCS cross-sections). Lattice QCD analyses are based on the 
7 —{ \ operator product expansion and lead to information on the lowest few Mellin moments of the GPDs which correspond 
\& • to non- forward matrix elements of twist- two operators. Again most studies focus on the proton [HI, O Hi ' but 
the GPDs of other hadrons are equally accessible in the lattice approach (unlike in experiment). In particular, recent 
studies by the QCDSF collaboration [15[ have investigated the GPDs of the pion. 

Here we investigate the vector and tensor GPDs of the pions from the point of view of chiral perturbation theory 
applied to lattice QCD, studying the quark mass and lattice volume dependence of the meson matrix elements of 
J^" 1, twist-two operators. Related studies of the vector GPD pion matrix elements have been presented in Ref. [l6|, [l?], [H, 
HHMm, however we extend those results to partially-quenched chiral perturbation theory (appropriate for lattice 
calculations with differing sea- and valence- quark masses). We also consider the tensor twist-two operators (also 
treated recently in Ref. [22|]) and include the effects of the finite volume to which lattice simulations are necessarily 
restricted (the Lorentz non-invariance of the lattice boundary conditions introduces some novel issues in the resulting 
extrapolation that we highlight). The finite volume calculations in partially-quenched chiral perturbation theory are 
particularly relevant to the ongoing lattice calculations reported in Ref. [l5| (see [Hf for earlier work in the forward 
limit). 

In the following section we provide our notation and conventions for the pion GPDs and twist-two matrix elements. 
In Section IlIII we turn to the effective field theory description of these objects before presenting our infinite volume 
results in Section IIV1 their finite volume analogues in Section [V] and a concluding discussion in Section IVII Various 
aspects of the finite volume forms are relegated to the Appendix. 

II. GENERALISED PARTON DISTRIBUTION OF THE PION 

A. Pion GPDs 

The generalised parton distributions of the pions (here we restrict our discussion to SU(2)) are defined by matrix 
elements of light-cone separated bi-local currents. Specifically, the vector GPDs H A , and the tensor GPDs, E£, are 
given by 

(^(P')l?(-fn)^7-^(|n)|7r j '(p)) = / dye*"** H A (y,£,t)u ■ Ptr [rW] , (1) 
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and 

(ttV)!^ (~|n) tHu^u^ (~n) ^'(p)} = y dye- iyzu - p E A {y^,t)u ■ Pr ■ Atr [rVV] , (2) 

respectively (for simplicity, we have suppressed the gauge links that render these matrix elements gauge invariant). 
Here the average four momentum of the incoming and outgoing pion states is P = h(p + p') and the momentum 
transfer is A = p' — p, u is a light-like vector (it 2 = 0) and r is transverse (r • u = r ■ P = 0). The GPDs are functions 
of the three variables y, £ = ^^j^ an d t = A 2 . Finally, ta — (1,Ta) for QCD. 



B. Twist-two operators 



As is evident from the forms of the GPDs, there are two towers of local twist-two quark operators that have 
non- vanishing pion matrix elements. These are given by 

Of = ... u^O'X -^ =l\>u . 7 (iu- D) n TAi> , (3) 

and 

Ofy = u a rp Ufll . . . u^O^X^ = ^ a al3 u a r {i u ■ D) n r A ^ . (4) 

The operators O^ 1 " and Oj.^ are of fixed twist (= dimension - spin) and transform irreducibly under the 
Lorentz group. Matrix elements of the vector twist-two operators in Eq. ([3]) give moments of the quark distribution 
in the pion in the forward limit. There are also two towers of purely gluonic operators at twist-two (see for example, 
Ref. [10|]) that have non-zero pionic matrix elements. For the purposes of our current discussion we note that the 

vector case has the same transformation properties as O ™^ above while the tensor case is beyond the scope of this 
work. 

These operators transform as (3, 1)0(1, 3) (isovector, A = 1, 2, 3) or (1, 1) (isoscalar, A — 0) under SU(2)l xSU(2),r 
rotations. The tensor twist-two operators also have non-zero matrix elements but vanish in the forward limit and 
belong to the (2,2) © (2,2) representation of SU(2) L xSU(2) fl irrespective of the flavour index A. In the SU(4|2) 
partially-quenched QCD case (where additional valence and ghost quarks are introduced), these operators are extended 
by the replacement of t a by f n , f„ and tt in the isoscalar- vector, isovector- vector and tensor cases. These matrices 
are somewhat arbitrary [2J, |25|, [26| , but for dcfiniteness we choose: 

f = diag(l, 1, 1, 1, 1, 1) , f 3 = dia,g(l,-l,qj,q k ,qi,qj+qk-qi), f T = diag(l, y, q 3 , q 4 , q 5 , q 3 + q 4 - q 5 ) , (5) 

which reduce to the usual Pauli matrices in the QCD limit and transform in the corresponding representations of the 
enlarged flavour group. 

The local twist- two QCD operators are simply related to those in Eqs. |T]) and @ and it follows that 

{n\p')\O^W{p)) = H A +1 (^,t) (u-PT +1 tr[fW] , (6) 
<7rV)l4 n ik>)> = *T in+ ifot) (r-A) (u ■ P) n+1 tr[fW] , (7) 

where H A +1 (£, t) = £ dyy n H A (y, £, t) and E A n+1 (£, f) = j\ dyy n E A (y, & t). 

Discrete symmetries and the approximate isospin symmetry of QCD constrain the pion GPDs. Time reversal 
invariance demands H A (y,^,t) — H A (y, —£,t) and E A (y,£,t) = E A (y,—£,t). Under charge conjugation (C), both 
the vector and tensor operator transform as [13, [27j : 

CO^C- 1 = (-l) n+1 0%\ (8) 
CO^C- 1 = {-l) n+1 { ^ A . (9) 

Using this, it can be shown that the isoscalar {A = 0) vector matrix elements vanish for even index n = 2k and the 
isovector (A = 3) vector matrix elements vanish for odd n (additional complications arise in the SU(4|2) case). For 

odd index, the isoscalar- vector matrix elements are parameterised in terms of generalised form factors A ® ^ and C„°^j 



as: 



(^(P')\O^W{P)) = 2sJ k J2(u-Af(u-P) 2k ~ 2 U^Jt) + (u-Af k C^Ht)) , (10) 
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while for even index, the isovector-vector matrix element can be parameterised as: 

<^(P')|ofV0P)> = 2ie« fc £( W -Af (u-P) 2k - 2l+1 A^ +li2l {t). (11) 

1=0 

Similarly, parameterisations of the tensor operator matrix elements are given by 

k 

<7rV)|0g + V(P)> = 2S l ^^(u-P)(r-A)(u-Af(u-P) 2k - 2l+1 B^ k % l (t) 1 (12) 

i=o x 

^(p')\O^W[p)) = 2 l£ ^^-l( U .p)(r.A)( U .Af (u-P) 2k - 2l B^l 2l (t), (13) 

I=i x 

for both even and odd index. 

The generalised form factors, and B^'^ A \ are related to the generalised parton distributions in Eq. (fTJ) and 
© as 



/l n-l 
-1 ■ n 

J— U, even 

/I n 
dyy n H (3 Hy,U) = E ("^ A nli J (*) . 
-1 o— n 



j =0, even 



and 



/I n 
"1 n — n ouot, 



(14) 
(15) 

(16) 



j— 0, even 



respectively. C^^t) is the D-term form-factor (28j . 

III. EFFECTIVE FIELD THEORY 

The small t behaviour of the GPD form- factors can be reliably described by the low energy effective theory of QCD, 
chiral perturbation theory (x?T) [2^, [3(|. The extensions of this to partially-quenched QCD 1 (in which valence 
and sea quarks have different masses as appropriate for many current lattice calculations), partially-quenched xPT 
(PQxPT) [H, HH, is well known and here we simply highlight the relevant pieces of the Lagrangian and discuss the 
operators that contribute to the twist-two matrix elements. We primarily focus on a partially-quenched theory of 
valence (u, d), sea (J, I) and ghost (u, d) quarks with masses contained in the matrix 

rriQ = diag(m„, m dl rrij,mi,mu, mj) , (17) 

where m- j = m Uj d such that the path-integral determinants arising from the valence and ghost quark sectors exactly 
cancel. 

A. Lagrangian 

At leading order the PQ^PT Lagrangian is given by 



f 2 f 2 r i 

= y str [D^Et^S] + A^jstr m Q ^+m^ Q Y. , (18) 



1 We do not discuss quenched QCD in which sea quarks are omitted as it has no connection to physical observables except in the large 
N c limit [H. 
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where the pseudo-Goldstone mesons are embedded non-linearly in the coset field 



E = exp 



/ 



(under a chiral rotation, E — > LURJ) with the matrix $ given by 

4> 



and 



M x t\ 



(19) 



(20) 



M = 



(f]u 
7T~ 



""J 



IT 

Vd <Pdj <Pdl 

\4>iu 4>id 4>ij m 



m = 



7T ' 
Vd 



b dd 



hi 



(21) 



The upper left 2x2 block of M corresponds to the usual valence-valence mesons, the lower right to sea-sea mesons and 
the remaining entries of M to valence-sea mesons. Mesons in M are composed of ghost quarks and ghost anti-quarks 
and are thus bosonic. Mesons in x contain ghost-valence or ghost-sea quark-anti-quark pairs and are fermionic. In 
terms of the quark masses, the tree-level meson masses are given by 



Ml. =M%. 



= A 



K)i, + (™q) 



(22) 



where Q ~ (u, d, j, I, u, d). The decay constant is normalised as / ~ 132 MeV. Additional terms involving the flavour 
singlet field, str[<I>] are not relevant here; in both PQ%PT and aTT the singlet meson acquires a large mass through 
the strong U(1)a anomaly and can be inte grat ed out, leading to a modified flavour neutral propagator that contains 
both single and double pole structures [3J, [35| . 



B. Twist-two operators 



Twist-two operators have been studied quite extensively in various low energy effective theories. A number of studies 



19 ,l2C,l21l b ut the relevant operators also 



(these studies have also been 



have focused on pionic matrix elements of twist-two operators [16, 
contribute in numerous studies of nucleon matrix elements [22], 1261127 
extended to the nuclear setting in Refs. 

We first focus on the vector operators, Eq. ([3]). To perform the matching of these operators to those in x?T it is 
useful to make the separation: 



O 



(m) _ Q(m) 



A,L 



O 



(to) 
A,R > 



(23) 



such that O 



(m) 
A,H 



Qht^(u ■ 7)(w • iD) m 1 qn for H — L,R, where ql,r = [(1 T 7s)/2](7 are the left (right)-handed 
quark fields which transform as qL — > Lq^ and qn — » Rq^ under the action of SU(4|2)lxSU(4|2)^. To construct the 
EFT operators, it is useful to treat f£ R as a spurion field that transform under global chiral rotations as: 



f^-Lf/Lt, T&^Rf&tf (24) 

(the spurion fields take a vacuum expectation value of ?a). This promotion renders the QCD operators, OT 1 ', 
invariant under chiral rotations. For these operators to have the correct charge conjugation properties, Eq. (|12p . 



requires f^ R —* (t£ r) T since E A E T and m q A (m,) T '. 

At leading order, the EFT operators consistent with these transformation properties are constructed from a single 
insertion of R and arbitrary numbers of E/E^ pairs. The first two terms in this tower are given by: 



m+l 



°a!i - E (<* str[^(D fc S)(D m - fc + 1 Et ) ] + a ^ s tr[? A (a k ^)(a m - k+1 X)} 

fc=0 

jl+j2+j3<m+l 

E 



{ Tl = V a A ' 2 



A,2 "mjlj2j3 
h,h ,33=0 



str[7U (a Ji E) (n h ) (□« S ) (□ 



-iTO+l-Ji -32—33 vt 



E t ) + (E^E t )] 



(25) 
(26) 
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where □ fc = {iu ■ d) k . Charge conjugation requires that a. 



A, la 



A,l i 
a rn,k(- 



-l) m+1 , and parity, E — > ET, implies that 



A,l 



,(m) 



(-1) 



m+1 



A,l 



significantly simplifying O^j . The second type of operator, , is not independent of 



(m) 



^ at infinite volume. If there are non-zero numbers of derivatives on three or four coset fields (E*-^), this operator 
may contribute through diagram (b) in Figure [1] below, but will be proportional to powers of u ■ k where k is the 
integration momentum. For even powers this will produce overall factors of u 2 which vanish; for odd powers, the 
integrand will be odd and hence vanish upon integration. Consequently, these operators only contribute in FigJTJb) 
when all derivatives act on the external pion fields. However using integration by parts and eliminating operators 
with derivatives on more than two meson fields such terms can be rewritten in terms of Operators involving 

six or more coset fields can similarly be eliminated. 

Thus the final form for the vector twist-two operators we use is: 



O 



(m) 



f2 m 



4 ^m+ijy 

j—0, even 



(-□) J '{str f A T,n m+1 - j tf +str fU E t 8 ro ~ J ' +1 E } 



(27) 



where we find it convenient to express the result in terms of the forward-backward derivative □ = □ — O = iu - (d — d) 
and time-reversal invariance limits the sum to even values. 2 At finite volume some of the operators we have neglected 
will contribute as Lorentz symmetry is no longer preserved. Here we ignore such terms but they result in additional 
complications in the extrapolation needed for lattice data as discussed in Section fVl 

Construction of the tensor operators is similar to that of the vector operators, however in analysing the trans- 
formation properties of the QCD operators it is necessary to introduce additional spurion fields, ff R and f RL , that 



transform as t£ r 



Lff R RK 



and t£ l 



Rt£ r L' such that 



,j-i(m) - -A ■ 



{iu ■ D) m q R + q R fi L tu^a^(iu ■ D) m q L , 



(28) 



is invariant under chiral rotations. 

As in the vector case, there is a tower of operators at leading order in the EFT with arbitrary odd numbers of 
E and E^ fields and a single insertion of ff R RL consistent with the requisite transformation properties. A similar 
discussion to that above for the vector operators greatly simplifies the structure of these operators and shows that 
only the operators involving three coset fields contribute to the single-particle matrix elements we are considering at 
next-to-leading order in infinite volume. Using charge conjugation and noting the QCD operator is antisymmetric 
under the interchange u <-> r, this relevant operator can be written as: 



°T\A = y E Wi J (-nr{(-Q)str[f T {E(Et3™+ 1 ^E) + (E0 m + 1 ^St)E} 



(29) 



j—0. even 



(-□)str T T {E (EtBB" 1 -^) + ("e E f ") E} +(E^E t )} 



where □ = (ir ■ d). 

Unlike the vector operators, the tensor operators involve a single set of low energy constants (LECs) because they 
belong to the same chiral representation regardless of the choice of flavour structure. In the SU(2) case, the super- 
traces in the operators of Eq. (|27[) arid Eq. (|29[) reduce to ordinary flavour traces and the various matrices are now 
2x2, but the form of the operators and the LECs that appear are otherwise unchanged. 

At the next-to-leading order (NLO), 0(p 2 ), vector and tensor operators which will contribute at tree-level are 
generated by combining the leading order operators above with insertions of d 2 , d a . . . d a or by substitution of the 
quark mass matrix mg for a coset field, E. The explicit forms of these operators are not required here however they 
generate polynomial dependence on the quark masses and t. 



IV. INFINITE VOLUME RESULTS 



At leading order, the moments of the pion GPDs receive tree-level contributions from the operators in Eqs. (|2T|) and 
(f2T))) above. At next-to-leading order contributions come from tree-level insertions of the 0(p 2 ) operators discussed 



2 For notational convenience, we use the same symbol to denote both the underlying QCD operator and that in the effective theory. Note, 
however, the EFT operators match not only the leading twist QCD operators but also higher twist operators of the same quantum 
numbers and also to purely gluonic operators in the isoscalar cases |42H . 



6 




FIG. 1: One loop diagrams contributing to the mesonic matrix elements of twist-two operators. The filled circle denotes a 
vertex from the chiral Lagrangian while the crossed circle corresponds to the twist-two operator. Diagram (c) denotes the wave 
function rcnormalisation. 



above and from the one-loop diagrams involving the leading order operators shown in Fig. [TJ The higher order 
operators lead to polynomial dependence of the GPD form-factors on m q ~ M 2 (M is a Goldstone meson mass) and 
t whilst the loops generate non- analytic dependence on these quantities. Therefore at NLO, the vector GPD form 
factors will have the form: 

4flu« = + A ^+h M2 + 4f+iV + loop contributions, (30) 

where A^ +l ^ is the bare matrix element determined in terms of the leading order LECs a^ +l j , and ■ and 

v4m+i . can be expressed in terms of linear combinations of the LECs accompanying the various NLO operators and 
absorb divergences from the loop contributions (in general this renders these terms renormalisation scale dependent). 
A similar expression holds for the tensor GPD form factors. 

In the following subsections, we present results for the various matrix elements in the SU(4|2) isospin limit case, 
jTid = 7Ti u and to; = m j. In this limit, only the valence- valence meson mass, M uu = M„, and the valence-sea meson 
mass, M 2 j = M 2 — ^S 2 j, enter. Here S 2 j/2X is the mass splitting between u and j quarks and results in SU(2) are 
easily obtained by setting 5 U j — ^ (these results are given in Appendix |B|) . 



A. Isoscalar-vector operators 

In the partially-quenched theory, the pion matrix elements of the isoscalar-vector operators have the form 

m 

(n+\0% n) \n+) - 2(l-(-lD E «i°Vi J [(2^-< l ^ +1 (A-^-(A- M r +1 ] 

j— 0,cvcn 



(31) 



(l-(-l) m ) 

48/ 2 7T 2 

+3 (l-a m+1 - j )lo, 



A) m+1 J2 



a 



(0) 



' fl2(l-o^H)^ 



,3=0, even 

M 2 + \t(a 2 -l) 



Mi + lo 



4Af2 + t {a 2 - 1) 

M 2 ul + \t(a 2 -l) 



x (t ({j - to - A)a 2 - j + to + 8) a m+1 - j + 4 ((j - to - 2)a m+1 - j + l) M 2 j + t (3a 2 - 7)) | , 

for to odd and vanish for to even. Note that above and in what follows, we have suppressed analytic dependence of 
matrix elements on M 2 , M 2 ^ and t. The first term in the braces arises from a double-pole propagator in Figure QJ a) 
and gives enhanced quark mass dependence. The QCD limit is easily obtained by setting S U j — and M U j — ► M n . 
Decomposing this result leads to the following structure for the GPD form factors (to odd): 



4 (0) (t\ - A [ "> 4- A [v ' a > M 2 -i- A ^ u ' a> M 2 4- A ylJ '"> t for n < m 
A m+l,j K 1 ) — A m+l,j + A m+X,j 1V1 Tt + A m+l,j 1VI uj + A m+l,j t -> J - m i 



|(0) 



,(0,o) 



(0,o) 



I (0,6) 



(32) 
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with no non-analytic dependence, and (m odd) 



r (o) ,1 



"'- 1 ' m+1 96^ 2 / 2 7-i "~ 1 4Af 2 + t (o2 _ i) 
+((3a 2 -7)t + 4M 2 i )lo" " ^ 



3M 2 log 



Mh - 4(1 - a 2 ) 



1 



1927T 2 / 2 



j= 0,cvcn 



■x f 12^-M 2 
[ 4M 2 +t(a 2 - 1) 

f 



3M 2 log 



(33) 



(1-a 2 ) 



+ i(3a 2 - 7) + 4Af 2 - + 4(m - j + 1) - -(1 - a 2 ) lo. 



A-f^.-Kl-a 2 ) 



where the LECs and the bare form-factors are related by C^' +1 = — 



(0) 



j=0,even m+1 J 



and A (0) - 2 m - j+1 a ( - 0) 



The ir~ and 7r° results are identical and, with appropriate changes in normalisation, the QCD limits of these results 
agree with those in Refs. [H, Hl[ using integration by parts and noting that J2j even ^ '^m+i j = ( tne resu lts for the 
gravitational form- factors, A2,o(i) and C2(i), also agree with previous calculations [43|). In the partially quenched 
case, this form-factor leads to an enhanced divergence in the pion gravitational radius ~ S^/M 2 as opposed to lnM^ 
in QCD. 



B. Isovector- vector operators 

The pion matrix elements of the isovector- vector operators have the form 



( 7 r+|Of ) |7r+) = (l + (-lD £ 2(2P.<^+ 1 (A.^(l-^M 2 J log 



j =0, even k 



M 2 ■ 



da a m - 1+2 log 



M^.-|(l- a 2 ) 



2M 2 .log 



M 2 



(34) 



for even m. This leads to the following structure for the GPD form factors: 



A {3) I 1 
m+lj I 8tt 2 / 2 



Ml- 
" 3 log 



M 2 



A^l,^ + + A^t, for j < m - 2 , (35) 



r(3,o 



(3,6) 



and 



^m+l,m(*) 



4 (3) ! I 



M ui I M u, 

— "Mog I — ^ 



4(3,0) M 2 , J(3,o) M 2 , A\0,U) , 



E 

j =0, even 



y- m A 



(3) 

771+ l,j 



8tt 2 / 2 



M 2 log 



.(3,6) 



M 2 ■ 



da a m - j+2 log 



AT 2 - - 4 



(l-« 2 ) 



(36) 



(3) 

using a> m ' +hj 



23- 



lj ^m+ij- The 7r results are related to these by factors of —1 and those in the 7r° vanish. 

The LECs A^ a) and i^ a) vanish by current conservation and A^l = 1. These results can be shown to agree with 

Ref. [H, [2l| (and earlier results in the case of the vector- isovector form- factor, A^l(t)) using integration by parts and 
noting the different normalisations. 

In the partially-quenched theory, isospin is not a good quantum number (the SU(4|2) adjoint matrices are given in 
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Eq. and the odd-m matrix elements are also non-zero. These take the form (tt and ir° are identical to the 7r + ) 



,(3) 



< 7 r + |Or ) k + ) m odd = {± ^jp {qj+qk ) J2 « m+llj 

j— 0, even 



(37) 



(it • A) r 



24 



da 



8 (2a> + 3 (a 2 - l) a m+1 ) A/ 2 (5+ 2 4 (a* - a m+1 ) M%5, 



5 2 . 



(4M 2 + <(a 2 - l)Y 



AM 2 + t (a 2 - 1) 



/A/ 2 + ±i(a 2 -l)\ , x /A/2, + k(a 2 -l)\ , 
+ log [ * 4 J ^ ) G(A/ 7r ) - log ( u ° y i G(M UJ -) 



+ [(« • A) m+1 (2P ■ ur^(u ■ AY] (fog Ml - log i^k j M 2. j I , 

where G(A/) = (i (3a 2 - 7) a J ' + t ((j - m - A)a 2 -j + m + 8) a m+1 + 4 (ot j + (j - m - 2)a m+1 ) A/ 2 ) . Note that 
this matrix element vanishes in the QCD limit and the sea-isospin limit where qj = —qu, in which foreseeable lattice 
calculations would be performed. Consequently, we do not present the form-factors that result. 



C. Tensor operators 

For m odd, the matrix elements of the tensor operator are given by: 

m ( 

<7r+|4 m) |7r + ) = 8(l-(-l) m )P-wA b m+1 , j (2P-u) m - j (u-Ay\(l+y) 

1 



48/ 2 7T 2 



j*=0,even 

2(3(g 3 + 94) + 5(l + t/))log 



— ^- log ' — ^ 



12tt 2 / 2 



H 2 j 113 



(38) 



Mi 



3 log -i ((2(g 3 + 94) - (1 + y))M 2 + (1 + y)^.) 



which simplifies considerably in the QCD limit. These matrix elements vanish if the charge matrix is isovector in 
both the sea and valence sectors. Matrix elements in the ir~ and tt° are identical. 
For m even, the matrix elements are given by 



(7r+|0^ m) |7r + > - 2(1 - y)P ■ ur ■ A ]T b m+1 A - 4(u ■ A)-? ((p • u)™-' - 2(2P • u) m ~ J ' + ((p + A) • u) m ~ j ) 

j— 0,evcn I, 



A/ 2 - /A/ 2 - 
— log ' 



12tt 2 / 2 



48/ 2 7T 2 



Slog ( ^f) (A/ 2 + _ io log ^ ] A/ 2 



i,/,Ai ' ' da(4A^+t(a 2 -l)) 



16/ 2 7T 2 J j 



1 - a 



m— j 



2a r 



1 + a 



m— j 



log 



A/ 2 3 + ^(a 2 -l) 
M 2 



(39) 



vanishing in the isoscalar case. The tt matrix elements differ by an overall sign and the 7r° matrix elements vanish. 
The results in Eqs. ([38|) and (139)) are easily converted into the various GPD form factors. For m odd, 



R T ft) - R J 'W 



_ R T,(0) J (i+y) 



f 



g3 + g4 A/ 2 log 
16tt 2 / 2 



— V±- log 

8tt 2 / 2 S ^ fi 2 

Ml " 



M. 



7T ^UJ 

16tt 2 / 2 



93 + 94 „ 2 



16tt 2 / 



A/ 2 ■ log 



A/ 



2 f 2 «J 



B 



T,(o) 



M'i + bL^Km 2 ^ +b. 



t. 



m-\-l,j 7r ' m+lj uj ' m+ljj 5 



(40) 



for all j. Here B^'j ^ = 8 2 m ^ A x b m+1:j (m odd). 
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While for m even, 



s m+i,i(*) - B 



T,(0) 1 - V 



m-\-l,j 



for j < m, and 



167T 2 / 2 



Mi 



log 



- 2M 2 ■ log 



R r,(a) M 2 



f>T,(a) M 2 



r>T,(b) 
JD m+l,j i i 



B 



T,(0) 
m+l,m 



i 



y 



16tt 2 / 2 



T,(0) 



(M^ + Slj) log 



Ml 



2M 2 fog 



M 2 - 



j— 0, even 



da M, 



t / 2 



log 



M 2 - 



-I? 



T,(a) 



Mf+B, 



T,(a) 



,M„. 



T,(6) 



where = -8A X 



m+l,m 7T ~t~ m+l,m uj ' m+l,m ' 
Efclu,even ( 7- fc ) 2 ~ J Wl.fc ~ 2 m ~ J 6, 



771+ 1, J 



(41) 



(42) 



(m even) is the leading order result for the form- 
factor. In the QCD limit, these result reproduce those of Ref. [22j [Eqs. (106) and (128) therein]. As these results are 
insensitive to the sea-quark charges, these form-factors can be calculated in lattice QCD without disconnected quark 
loops. 



V. FINITE VOLUME EFFECTS 



The space-time lattices used in numerical simulations are of finite extent by necessity, and consequently, lattice 
results differ from those of the physical (infinite volume) world even at the physical quark masses. 3 For sufficiently 
large volumes, these effects can be incorporated into the effective field theory approach, allowing infinite volume results 
to be extracted from the finite volume (FV) lattice simulations. Here we shall consider a hyper-rectangular box of 
dimensions L 3 x T with T > L as appropriate for most current lattice calculations. Imposing periodic boundary 
conditions on mesonic fields leads to quantised momenta k = (ko,k), k — = ^ (ji,j2, J3) with 6 Z, but ko 
treated as continuous. On such a finite volume, spatial momentum integrals are replaced by sums over the available 
momentum modes. This leads to modifications of the infinite volume results presented in the previous section; the 
various functions arising from loop integrals are replaced by their FV counterparts. In a system where M^L 3> 1, 
finite volume effects are predominantly from Goldstone mesons propagating to large distances where they are sensitive 
to boundary conditions and can even "wrap around the world" . Since the lowest momentum mode of the Goldstone 
propagator is ~ exp(— M V L) in position space, finite volume effects will behave as a polynomial in 1/L times this 
exponential provided no multi-particle thresholds are reached, that is for t < 4Af 2 . For t > 4M 2 , volume effects 
that are polynomial in inverse powers of L are expected, however for realistic lattice calculations, such momentum 
transfers are too large to be described in xPT and we neglect the resulting complications in our analysis. 

The finite volume effects in the diagrams of Fig. [ljb) and (c) arising from the operators that contribute at infinite 
volume are well known. However the effects in Fig. [IJa) are made more complicated by the non-trivial numerator 
structure of the integral/sum and the requisite details for their evaluation are given in Appendix [X] The final forms 
for the finite volume versions of our results are given below. In most cases, these replacements of spatial momentum 
integrals by sums would complete the calculation of finite volume effects. However for the GPD form-factors, there 
are additional complications that arise from operators whose contributions vanish in infinite volume but which give 
non-zero contributions at finite volume. To see this, we concentrate on the vector case and reconsider the operator in 
Eq. (|26p (note there is some redundancy of terms in this operator). Naively, this operator would contribute through 
diagram (b) in Fig.[T] As discussed in Section [III[ there is no contribution at infinite volume as the each of the pion 
fields coming from this operator must have a non-zero number of derivatives acting on it. Combining u 2 = (or 
equivalently, the definite Lorentz transformation properties of the twist-two, spin-n, operators) and the symmetric 
range of integration, terms with derivatives inside the loop in Fig. [ljb) give no contribution. 



3 Lattice artifacts from the discretisation of space-time also influence lattice results. We do not discuss these here and assume a continuum 
extrapolation has been performed. 

4 In principle, finite volume effects can also be computed for M^L ~ 1 but A X L 2> 1 where mesonic zero-modes become enhanced 
[44L |45| . l4al . These calculations are beyond the scope of this work. 
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In the finite volumes discussed here, Lorentz symmetry (or 0(4) symmetry in Euclidean space as relevant in lattice 
calculations) is explicitly broken by the imposition of the boundary conditions (in a L 3 x T periodic box, arbitrary 
Lorentz boosts do not leave the system invariant). As a result, the arguments used to discard the operators in Eq. (|26|) 
break down (-^ J dqo Y]^ r^z^ru 7^ at finite volume) and their contributions, which must vanish as L — > 00, need to 
be incorporated to correctly describe finite volume lattice calculations Consequently, one must include the functional 
dependence arising from operators such as those in Eq. (|26p in any fit to lattice data and determine the relevant 
combinations of the a^t,' 2 ■ ■ as well as the a^t,^ „ before discarding the former (we have no interest in these LECs 
for the matrix elements we want to extract, however they would contribute in more complicated matrix elements such 
as two pion matrix elements of twist-two operators) in extracting the infinite volume result. In general this is a very 
intricate task as the functional dependence produced by these operators in the contribution of Fig. QJb) to the matrix 
element of OT^ is 



E 

1—2. even 



1 x k^ 1 k^ 1 — tr 
7n,I«/ii • • I dk jjJ2 k 2 '_ m 2 ( 43 ) 



where the parameters j n j are linear combinations of those that enter in Eq. (|26p . 5 However for the low moments 
(to < 4) that are accessible in current lattice calculations, these subtractions can be performed. 

A number of additional aspects of Lorentz symmetry violation at finite volume are worth noting. Firstly, these 
effects did not contribute in the nucleon twist-two matrix elements at finite volume at NLO [2(|, but will contribute 
at higher orders in the chiral expansion in a similar way as they enter here. We also note that moments of distribution 
amplitudes (given by meson to vacuum matrix elements of related quark-bilinear operators) will suffer from similar 
complications at finite volume and the absence of non-analytic quark mass dependence at NLO found in Ref. [2(| 
at infinite volume will not persist. Secondly, as shown in Ref. [45j |. no new operators can appear at finite volume 
otherwise their LECs would necessarily depend on L~ l , an infrared scale. In the case described above, the operator 
is present both at finite volume and infinite volume, but doesn't contribute to the matrix elements we consider in the 
latter case. A final related effect of the finite volume is that Lorentz symmetry can no longer be used to decompose 
the matrix elements of the twist-two operators. At finite volume, additional structures that violate Lorentz symmetry 
can appear in the decomposition of the matrix elements in Eqs. (fT0f - (fT3)) 6 . The form factors of such terms vanish as 
L — > 00 and we ignore them here but care must be taken to correctly extract the form- factors in Eqs. (fTUHT5)) without 
pollution from these additional Lorentz non-invariant form-factors. 

A. Finite volume matrix elements 

With the preceding remarks in mind, the finite volume versions of the results of the previous section are given by 
the following matrix elements: 

m ( 

(7r+|^ m) | 7 r+)?V odd = 2(1 -(-I)™) J2 ^ +1 A((2P.ur +1 -^u-Ay-(u-Ar +1 ) (44) 

j— 0,evcn ^ 



(u-A) 



m+l 



- 2 (2Ml - S 2 U j + 2b) Jr 7 , 0i0 , 



3/ 2 

-\-3M^ (^0,0,0,0,0 + (-^0,0,0,0,0 + ^o,o,o,o : o)) + 4 (^0,0,0,1,0 + ^0,1, 0,0,0 
m+l— j 

£ (-+ 1 f-')2*(u.A)-*(-2(2JI^-^ + 2t)A2j, i0j0i0 



fc=0 



+3M 2 (-^0,0,0,0 + $uj (^™oao,o + ^"ao,o,o)) + 4 (^,0,0,1,0 + ^M, 0,0,0) ) 



5 Care must to be taken to define a suitable regularisation prescription. 

6 In the EFT analysis of the matrix elements under consideration, such terms do not appear until higher orders (two-loops) in the chiral 
expansion, but even to the order we work, form-factors acquire dependence on \u\ 2 . 
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(Tr+lofV)™ = 2(1 + (-1D Yl ^ +1 A(2P-u) m +^(u-Ay(l-^ 

j — 0,even L 



2 

J2 "'">•' 



(45) 



2 m+l-j 
J k=0 



2 k (v ■ AY n+1 - k ( K uj + 2K U ° \ I 

z V u ^1 y^kfifi, 1,0 ^ z "fe, 0,1, 0,0^ f ' 



(7r+|Of l V + )Fv° dd 



(l-(-l)™) 
3/ 2 



(<&+?*) E ffl £+U U«-^r +1 [-( 2M '-^+ 2*) ^0,0,0,0 (46) 

j=0, even ^ 



(2Ml - S 2 uj + 2t) ^0,0,0,0 + 2 (M 2 + i) i^o,o,o,o 



+2 (M 2 + t) ^0,0,0,0,0 + SM^S^j (Kq^ + ^o,ofo,o,o + ^j^o,o 3 o,o,o) 

1 9 / _ k uu 4- K u i — K uu \ 

+ z ^0,0,0,1,0 A o,o,o,i,o + A o,i,o,o,o A o,i,o,o,oy 

m+l—j 
k=0 

+3M 2 <5 2 y (^S|o,o,o + ^feAo,o,o + ^«j^fc,o 3 o,o,o) 
+ z ^fc.o, 0,1,0 ~ A fe,o,o,i,o + A fe, 1,0, 0,0 A fe, 1,0,0,0 y 

+6 (2^'+ m+1 (P • • A)- 5 - (u ■ A) m+1 ) (J UJ - J„,„) 



2(1 - y)P ■ ur ■ A E 

.7=0, even 

16(1 -y) 



B r(o) 



'^( U -A)'(/'.«)"' »i, -2./,,, - /,,„) 



/ 2 



m 1 

E - r ■ A(u ■ Ay (iC+i-,,0,0,1,0 + 2^ +1 _ ji0 ,i,o,o) 



.7=0, even 



+(u-Ay+i (^,0,0,1,1 + 2^,0,1,0,1) 



m-j 



\m — k 



fc=0 



E( m ^)(2 fc+1 -1) («-A) 

<>• a (a • 2a; : :,.,„ ; ) 



^fc+i.o, 0,1,0 ' z "-fc+i,o,i,o,o y 



(47) 
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and 



<7r+|^ n) |7r + >FV dd 



8r-A J2 {2P-u) m+1 - j (u-Ayb m+ltj 

j =Q, even 



(V + 1) ( 1 - ^(2J«j + Ju.u - SljL UlU ) ) (48) 



^2 (93 + 94) ( Ju,j - Ju,u) 



16(«- A) ? 

3F - 



E WijE(v)2 1 («- a )" 



x < m • A 



j— 0, even 
' 2 5 2 - 



(2Af 2 - S 2 uj + 2t) 034^0,0,0,1 - 2(234^0,0,1,1 - 2034^,0,0,1 

+A/ 2 (Q34^o,o,o,i - 234^S,o,o,i - 3 (g 3 + 94) S 2 uj (^0,0,1 + ^fc,o 2 o,o,i) 
+3Q34'5 u;) ^fc i o 3 o ,0,1 ) — 2<2 34 (— i-K^fc,o,o,o,i + ^&,o,o, 1,1 + -^m, 0,0,1) 



•A 



(2M 2 - 5 uj + 2t) 034#fej 1 ,o,0,0,0 - 2O34iffe|x,0,0,l,0 ~ 2 034-?C+l, 1,0,0,0 

+M 2 . (034^+1, 0,0,0,0 — 034-^fe+i, 0,0,0,0 — 3(^3 + 94) S^j (-^fc+1,0,0,0,0 + ^fe+1,0, 0,0,0) 
+3034^^^+1,0,0,0,0) — 2Q 34 (— iiffe+i,o,o,o,o + Kk+i, 0,0,1,0 + ^fe+1,1,0,0,0) 



where Q34 = l + y + (73 + g 4 , Q 34 = l + y — (73 — q 4 and Q 3 ' 4 = l + y + 2q 3 + 2g 4 . Note that the entire last sum in the 
above expression vanishes identically at infinite volume. 
Here 



_ f dq 1 \^ 

Ja - b ~ / <>*■ r.a 2^ 



2tt L 3 ^ q 2 - Af 2 + ie ' 



L a ,b 



dq 1 \ - 
9^73 2^ ' 



2^ L 3 ^ (g 2 - M 2 b + ze) 2 



and 



where 



K 



fg _ ^/,ff;i,i 



i,a,b,c,d i,a,b,c.d> i.a.b,c,d i,a,6,c,d' 



A 



/91 _ /^/.s; 2 . 1 



i,a,b,c,d 



dq 1 \ - 



7^/32 r /,g;l,2 T^/t/3 r /,s;2,2 

i,a,b,c,d '^i,a,b,c,d' i,a,b,c,d ,Kj i,a,b,c,d' 

{u- q) l {q 2 ) a (p- qf{q- A) c {q-r) d 



2tt L 3 ^ (g 2 - M) g + ie)»i ((3 + A) 2 - M 2 + ie) 



(49) 



(50) 



(51) 



/s 



which in turn is simply represented in terms of the functions appearing in Appendix IA1 

Since the external momenta p = P — A/2, A, u and r occur in the various integrals/sums, it is not possible to 
decompose these matrix elements into form-factors without specifying to a particular m. 



VI. DISCUSSION 

The results for the quark mass and volume dependence of the moments of the pion GPDs presented above are 
useful in extrapolating lattice data at small t in the chiral regime (the range of volumes and pion masses for which 
the chiral expansion is convergent) to the physical world. Although subtleties not seen in simpler quantities arise 
at finite volume that complicate the extraction of the infinite volume results for the twist-two matrix elements, we 
have outlined a procedure by which the extrapolation can be performed reliably. The knowledge of the meson GPDs 
that can be obtained by combining our results with lattice calculations will be useful in comparison to, and as a 
complement to, experiment. 
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APPENDIX A: FINITE VOLUME FUNCTIONS 



The functions used in our results at finite volume can be built from the following basic structure 



Jw 2 (mi,m 2 , A;L) = -i j ^-j^^ 



2ir L 3 ^ (q 2 -m\ + ie)"i {{q + A) 2 - m 2 + ie) n ? 



1 f ST <9" 1+n2_1 

' L*r(m)r(n 2 ) J dx ^dk^ +n '- 1 

k 



k 2a k 11 ...k^ 



k Q + fc 



"1+112 



(Al) 



k =k 



where 

fc = \J \k\ 2 + 2(1 - x)k ■ A + (1 - .x)|A| 2 + xm\ + (1 - x)m\ , 
after performing the energy integral by contour integration. As an example, we find that 

Hi 



^^,0,0,0,1 — > u ...u r j k+111 u u ...u r <Jk,i,i 



k=0 



(A2) 



It is possible to write down a general expression for the multiple derivative in Eq. (jAip 7 but since we only require 
{ni,7i 2 } = {1, 1}, {1, 2}, {2, 1}, {2, 2}, it is easier to present each case explicitly. This leads to 



8L 3 



Ja,if(mi,m 2 ,A;L) = 

^"(rni^a, A;L) = 
where M = x(l — x)|A| 2 + xm\ + (1 — x^m^ and 



AL 3 J i-« 
(4a 2 - 8a + 3) f 1 



{l~x)A,M 
dx(l-x)lT- l J (l-x)A,M 
[l-x)A,M 



ia 3 - 36a 2 + 46a - 15) f 1 



16L 3 



dxx{l-x)I l l^ (l-x)A,M 



I^(z,m) = J2 



q" L . . . q' 



(\q + z\ 2 +m 2 ) 



2\P 



(A4) 
(A5) 
(A6) 
(A7) 

(A8) 



7 We may use the well known formula of Faa di Bruno 



dxn 9(m)-^-i^ m=ihl df[t)K in „ ) ' ( A 3) 



where K = £? =1 fc; 
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These sums can then be further simplified using the recurrence identity 



1 al /3-l 



z ln I l p'" ln 1 (z 1 m) 



2(0 - 1) dz in 

n 

H)"EE P ■ ■ ■ d * kztk+1 ■■■ zln ) W*. m) 



(A9) 



fc=o -p 



where V denotes a permutation of the orderings of the partial derivatives, d\ = ^j-, and the z"s, and the sum is over 
all such permutations. The remaining momentum sums have scalar summands and are simple to evaluate using the 
results of Ref. [ItJ and the tracelessness condition. Divergences can be regulated using Epstein-Hurwitz zeta-function 
techniques [4q |. 



APPENDIX B: QCD FORM-FACTORS 



For easy reference, the non-vanishing, infinite volume form-factors in the SU(2) theory are given here. 
In the vector-isoscalar case for to odd, 



and 



(Bl) 



m 

i: >p oj-myiW 



/^(O.a) , -=V(0,a)\ M 2 , /-r(0,6) , 



1927T 2 / 



(B2) 



j— 0,cvcn 



I daa ra - J+1 (i(3a 2 -7) + 7M 2 + 4(TO-j + l) ( Af 2 --(l-a 2 )) ) log 



M 2 -|(l-a 2 ) 
M 2 



In the vector-isovector case for to even, 



At\iAt) = A% +hj ( 1 - 7-^log ( ^ ) ) + + A mZi,j) M % + for J < to - 2 , (B3) 



I (3) 



M 2 



87T 2 / 2 



,(3,a 



(3,a) 



.(3,6) 



and 



4(3) m _ 4 (3) f i _ _^i np . f U , M (3 ' a) 4-4( 3 ' Q ) )M 2 +A (3 < b) t 



,(3) 



M 2 



m OJ — "i A \ 



8vr 2 / 2 

j ^0, even 

Finally in the isoscalar tensor case for m odd 



M>g 



tjT /.\ _ R T,(0) 



3M 2 



16tt 2 / 2 

while in the isovector tensor case for to even, 



B T (t) - B TA0) 



1 



272 l0 S 



16tt 2 / 



M 2 



M„ 2 m 
M 2 



M 2 



da a m ~ J ' +2 log 



M 2 -£(l-a 2 )i 



l°m+l.? + £> m+l,i) 1V1 -K n m+l,i l ' 



fj gT,(a) + £T,(«) )M 2 +B T,(&) 



for j < to, and 

nT /y> _ R T,(0) 

m+l,m V / m+l,rn 

l 

+ 16/ 2 7T 2 



16tt 2 / 



j =0, even 



/ jjT,(a) 6T,(o) \Af 2 ^R T '( h ) f 

(-"m+l.m ' - D m+l,m7 i "7r "T" - D m+l,m'' 





a" 


m— j 


) 


.2. 





log 



(B4) 



(B5) 



(B6) 



(B7) 
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